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Abstract 

In  this  paper  a  parallel  simulation  of  an  experimental  dual-mode  scramjet  combustor 
configuration  is  presented.  Turbulence  is  modeled  with  the  k-s  two-equation  turbulence  model  and 
a  7-species ,  8-equation  kinetics  model  is  used  to  model  hydrogen/air  combustion.  The 
conservation  form  of  the  Navier-Stokes  equations  with  finite -rate  chemistry  reactions  is  solved 
using  a  diagonal  implicit  finite-volume  method.  Using  about  3,120,000  grid  points,  the 
three-dimension  flow-fields  with  equivalence  ratio  4>=0.0  and  0.35  have  been  respectively 
simulated  on  the  parallel  computer  system,  obtaining  more  detailed  flow  properties  than  the 
experiment’s  results.  Wall  pressure  comparisons  between  CFD  and  experiment  show  fair 
agreement.  For  4>=0.35 ,  the  fuel-penetrating  height  of  the  seven  injectors  are  different  because  of 
the  effects  of  the  boundary  layer  and  the  shock  wave  in  the  combustor.  According  to  numerical 
results,  if  adjusting  the  locations  of  the  injectors,  the  combustion  efficiency  could  be  improved. 

Introduction 

Supersonic  combustion  ramjet  engine  (Scramjet),  in  which  the  combustion  process  can  be 
supersonic,  is  a  new  type  of  air-breathing  propulsion  device.  It  has  remarkable  perfonnance,  great 
potential  superiority,  and  a  good  prospect  in  the  field  of  military  and  civil  aviation  in  the  future. 
However,  the  experimental  study  of  scramjet  engine  puts  higher  demands  on  simulating 
capability  of  ground  test  facilities,  and  greatly  increases  the  test  cost  and  period.  Since  the  flow 
residence  time  within  a  combustor  is  very  short,  on  the  order  of  one  millisecond,  it  is  difficult  to 
measure  the  flow-fields ,  generally  only  get  the  wall  pressure  and  heat  flux.  Even  using  advanced 
measure  technology(Planar  Laser  Induced  Fluorescence),  which  is  expensive,  it’s  still  difficult  to 
make  three-dimensional  measurements,  hence  difficult  to  obtain  detailed  flow  properties.  With  the 
development  of  computer  technology  and  the  advancement  of  numerical  calculation, 
Computational  Fluid  Dynamics  (CFD)  can  provide  detailed  flow-field  characteristics  than 
experimental  results  and  has  emerged  as  an  extremely  valuable  and  cost-effective  engineering  tool 
in  combustor’s  design  and  analysis. 

Previous  work  in  the  area  of  combustor  analysis  and  predictive  capability  includes  the  early 
work  of  Billing1,  Billing  et  al2,  and  Waltrup  and  Billing3'4.  They  developed  correlations  based  on 
one-dimensional  and  two-dimensional  flow  analysis  which  included  such  drivers  as  incoming 
boundary-layer  thickness  and  maximum  combustor  pressure  rise.  These  correlations  perform  well 
on  simple  model  geometries  (axi-symmetric  or  two-dimensional)  but  are  not  generally  suitable  for 
more  complex  three-dimensional  problems  except  for  providing  directions  and  trends.  Recent 
computational  studies  of  the  three-dimensional  combustors  are  described  in  papers  by  Rodriguez 
et  al5  and  Riggins6.  The  computational  domains  include  the  use  of  jet-to-jet  symmetry  and  entire 
half-duct  modeling.  Calculations  in  each  case  are  typically  conducted  on  fixed  structured  grids  of 
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less  than  2,500,000  points  in  multi-zone  domains.  Among  all  the  numerical  simulations  the 
fuel-penetrating  heights  of  different  injectors  have  not  been  obtained. 

This  paper  is  intended  to  provide  a  parallel  simulation  of  an  experimental  dual-mode  scramjet 
combustor  configuration.  The  full  duct  for  the  isolator-combustor  is  then  modeled  in  order  to 
study  the  fuel-penetrating  height  of  the  seven  injectors.  The  grid  used  for  the  present  calculations 
is  211x41x121  for  the  isolator  and  211x81x121  for  the  combustor,  about  3,120,000  grid  points. 

Using  the  Navier-Stokes  equations  with  the  addition  of  conservation  of  chemical  species,  the 
three-dimension  flow-fields  with  equivalence  ratio  4>=0.0  and  0.35  have  been  respectively 
simulated  on  the  parallel  computer  system,  obtaining  more  detailed  flow  properties  than  the 
experiment’s  results.  Wall  pressure  comparisons  between  CFD  and  experiment  show  fair 
agreement. 

Governing  Equations  and  Numerical  Method 

In  general  curvilinear  coordinate,  the  full,  three-dimensional  Navier-Stokes  equations 
coupled  with  chemical  nonequilibrium  processes  in  a  nondimensionalized  conservation  form  are 
written  as 

8U  8E  dF  8G  if  dEv  8FV  dGv)  _ 

8t  8"t,  dr|  8C,  R 1  f  <9£,  dij  8C,  ^ 

where  U  are  the  conservation  variables,  E,  F  and  G  are  the  convective  fluxes,  Ev ,  Fv  and  Gv  are 
the  viscous  fluxes,  and  S  are  the  production  and  reduction  rates  for  the  nonequilibrium  processes. 

The  conservative  form  of  the  equations  is  solved  using  a  diagonal  implicit  finite-volume 
method,  which  approximately  solves  this  system  using  two  sweeps  of  a  point  Gauss-Seidel 
relaxation.  In  the  explicit  part,  the  inviscid  fluxes  are  computed  using  Steger- Wanning  scheme 
with  a  3rd-order  accurate  MUSCL  interpolation.  The  viscous  fluxes  are  evaluated  using  standard 
central  differences.  For  the  implicit  calculation,  the  code  uses  a  diagonal  algorithm,  which 
eliminates  the  expense  of  inverting  large  block  matrices  that  arise  in  chemically  reacting  flows. 

The  first  step: 

'1  +  A'(  Umax  +  Amax  +  Amax)  " 

C At-  \wu+ q  -  AtRHS 

The  second  step: 

1  +  AtU,,  max  +  rBmax  +  ACmaxf,)T/"  + 1  +  A/ 

i  AI+  l*f+ 1 +  Bj+ 1»”+ 1 +  A+  M  + 1)=[1  +  A<  A  max  +  A  max  +  Amax)]®  "  +  ‘ 

The  Chemical  Reacting  Model 

A  combustion  model  for  hydrogen/air'  mixtures  consisting  of  8  reactions  and  7  species, 
including  6  reacting  species  H2>  FFOx  O,  H,  OH,  O2  and  an  inert  specie  N2,  is  selected  for  our 
computations.  The  eight  reactions  are: 
l)H+02<^0H+0 
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2)  0+H2<=>0H+H 

3)  H2+0H-^H+H20 

4)  20H<=>0+H20 

5)  H2+M i <=>2H+M  i 

6)  H20+M2^0H+H+M2 

7)  0H+M3<^0+H+M3 

8)  02+M4<=>20+M4 

The  reaction  rates  are  taken  from  Evans  and  Schexnayder8. 

Boundary  Condition 

The  free-stream  is  supersonic  so  that  all  flow  variables  are  known.  Exit  conditions  are  set  as 
supersonic  extrapolation.  The  solid  wall  boundaries  are  modeled  as  no-slip,  adiabatic,  a  zero 
normal  gradient  of  pressure  and  fully  non-catalytic. 

Results  and  Analysis 
1.  Code  Validation 

Shock  induced  combustion  phenomena,  ranging  from  decoupled  to  coupled 
shock-deflagration  systems,  were  experimentally  investigated  in  the  mid  1960’s  and  early  1970’s. 
Calculations  of  this  type  of  reaction  are  very  demanding  in  terms  of  numerical  robustness  and 
accuracy,  since  the  reactions  usually  occur  very  fast  with  significant  energy  release  which  takes 
place  in  a  very  short  distance.  In  this  paper,  to  validate  the  code,  two  results  from  Lehr’s9 
experiment  are  reproduced  numerically.  The  two  cases  are 


Case 

Mixture 

Voo(m/s) 

Poo  (Pa) 

Too(K) 

1 

H2/Air 

2605 

42662 

250 

2 

H2/Air 

1685 

42662 

250 

All  cases  involve  a  sphere  having  a  diameter  of  15mm  moving  through  a  stoichiometric  mixture  at 
velocity,  pressure  and  temperature  indicated  at  the  above  table.  Calculations  are  perfonned  on  a 
45x80x41  grid. 


Experimentally,  Lehr  found  that  Case  1  resulted  in  a  coupled  shock-deflagration  system, 
whereas  Case  2  resulted  in  a  decoupled  shock-deflagration  system.  The  temperature  contours  of 
Case  1  along  with  the  experimental  shock  location  are  shown  in  figure  1.  Case  2  is  a  condition 
where  a  projectile  is  traveling  at  a  speed  lower  than  the  detonation  speed  of  H2/Air.  This  case 
resulted  in  a  decoupled  shock-deflagration  system  shown  in  figure  2.  These  two  calculations 
correctly  predict  the  location  and  the  shape  of  the  wave.  So  the  code  is  right  and  reliable. 

2.  Description  of  the  experiment 

This  section  briefly  describes  the  experimental  configuration  examined  in  this  investigation. 
The  scramjet  is  essentially  a  rectangular-section,  constant-width  duct.  On  both  upper  and  lower 
walls  there  is  a  backward  facing  step  and,  further  downstream,  a  3°  wall  expansion  only  on  the 
upper  wall.  These  two  features  divide  the  scramjet  into  three  sections:  isolator,  combustor  and 
expansion.  The  steps  have  a  height  of  3mm.  Injection  takes  place  95mm  downstream  of  the  steps. 
There  are  seven  injectors  on  the  bottom  wall  and  the  diameter  of  each  injector  is  1.2mm. 
Hydrogen  is  injected  through  the  orifices  at  sonic  conditions,  with  equivalence  ratio  0  =  0.35. 
Other  dimensions  can  be  obtained  from  figure  3.  For  convenience,  the  nominal  inlet  conditions  of 
the  vitiated  flow  are  provided  in  table  1.  Table  2  gives  the  injectant  conditions.  This  experiment 
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has  been  performed  at  CARDC. 

Table  1:  Nominal  inlet  conditions. 


M 

2.05 

Ph20 

0.258 

P 

0.341Mpa 

P  02 

0.21 

T 

1172K 

PiN2 

0.532 

U 

1417(m/s) 

Table  2:  Injectant  conditions 


M 

1.0 

Ph2 

1.0 

Pt 

2.3Mpa 

D 

1 ,2mm 

Tt 

300K 

0.35 

3.  Results  and  Analysis 

1)  0  =  0.0 


Without  injection  the  wall  pressures  at  the  bottom  on  the  center-plane  are  compared  in  figure 
4  with  the  data  of  CARDC ’s  experiment.  The  results  show  fair  agreement. 

2)0  =  0.35 

For  0  =  0.35 ,  figure  5  shows  the  wall  pressure  comparisons  at  the  bottom  on  the  center-plane. 
The  peak  pressure  is  almost  equal  and  the  results  show  also  good  agreement.  Figure  6  gives  the 
center-plane  mach  number  floods.  For  this  low-O  condition  the  whole  core  flow  is  supersonic. 
Figure  7  shows  the  wall  pressure  comparisons  of  the  numerical  simulations  between  0  =  0.0  and 
0  =  0.35.  For  0  =  0.35,  there  is  no  upstream  interaction  and  the  combustion  reactions  only  affect 
the  downstream  flow-field.  In  order  to  study  the  fuel-penetrating  height  of  the  seven  orifices,  the 
hydrogen-specie  flood,  the  temperature  flood  and  the  pressure  flood  of  the  cross-section  at  the 
center  location  of  the  orifices  are  shown  in  figure  8.  The  hydrogen-specie  flood  shows  the 
fuel-penetrating  heights  of  the  seven  injectors  are  different.  The  temperature  flood  indicates  that 
the  fuel-penetrating  height  of  the  two  side  orifices  is  low  for  the  impact  of  the  side  boundary-layer. 
At  the  same  time,  because  of  the  suppressing  of  the  shock  plane,  the  penetrating  heights  of  the 
three  middle  orifices  are  also  low.  According  to  numerical  results,  if  adjusting  the  locations  of  the 
injectors,  the  combustion  efficiency  could  be  improved 

Conclusions 

1.  Wall  pressure  comparisons  between  CFD  and  experiment  show  fair  agreement.  So  the  code  can 
be  used  for  studying  the  data  of  experiments  and  getting  more  detailed  properties  of  the 
experimental  flow  field. 

2.  Because  of  the  effects  of  the  boundary  layer  and  the  shock  wave  in  the  combustor,  for  4>=0.35 , 
the  fuel-penetrating  height  is  low  and  the  combustion  efficiency  is  low.  Adjusting  the  locations 
of  the  injectors  could  improve  the  combustion  efficiency. 

3.  For  0  =  0.35,  the  whole  core  flow  is  supersonic  without  upstream  interaction.  The  combustion 
reactions  only  affect  the  downstream  flow-field. 

Future  Work 

1 .  Unsteady  analysis  may  be  used  for  understanding  the  physics  of  fuel  air  mixing  and  for 
combustion  instabilities. 

2.  Consider  the  turbulence-chemistry  interactions. 

3.  More  complex  reaction  mechanisms  for  hydrogen/air  combustion  will  be  needed,  particularly 
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for  hydrocarbon  fuel. 

4.  Increase  the  convergence  rates  of  the  code. 
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Figure  1:  Coupled  shock-deflagration  Figure  2:  Decoupled  shock-deflagration 
system;  Temperature  contours  system;  Temperature  contours 
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Figure  3:  Outline  of  the  dual-mode  scramjet  experiment 
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Figure  4:  Wall  pressure 
comparisons  for  0  =  0.  0 
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Figure  5:  Wall  pressure 
comparisons  for  0  =  0.  35 
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Figure  6:  Mach  number  floods 
of  the  cnterplane  for  0  =  0.35 
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Figure  7: Wall  pressure  comparisons 
between  0=0.0  and  0=0.35 
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(a)  Hydrogen-specie  flood  of  the  cross-section  for  0=0.  35 
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(b)  Temperature  flood  of  the  cross-section  for  0  =  0.  35 
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(c)  Pressure  flood  of  the  cross-section  for  0  =  0.  35 


Figure  8:  Comparisons  of  the  fuel-penetrating  height  of  seven  orifices 
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